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National Aeronautics and Space Administration, Research Center, 

Hampton, Langley, Va., Grant NSG ll*tO, Computer Automation and 
Feedback System Design, June 30, 1976. 

SUMMARY 

Previous reports described a completely automatic computer program for the 
design of single input-output, single-loop feedback systems with parameter un- 
certainty, to satisfy time-domain bounds on the system response to step commands 
and disturbances. The inputs to the program are basically the specified time- 
domain response bounds, the form of the constrained plant transfer function and x 
the ranges of the uncertain parameters of the plant. The program output consists 
of the transfer functions of the two free compensation networks, in the form of 
the coefficients of the numerator and denominator polynomials, and the data 
on the prescribed bounds and the extremes actually obtained for the system 
response to commands and disturbances. This program was delivered and demon- 
strated to NASA personnel with reports describing the computer program. This 
report is therefore devoted to additional work since done. 


1. Explicit solution for optimum L(jto) . 

At some point in the program, there are prepared bounds il»(co) in the complex 
plane on the frequency response of the loop transmission function L(joo) -- see 
Fig. 1, i.e. L(jw) must lie outside the region described by » in order to 

satisfy the response specifications over the range of plant uncertainty. The 
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computer then seeks to find a rational function L(jco) which satisfies these 
bounds, using as a vital aid that the optimum L(jw) lies on its bound ip(w) 
at each ai . The technique used has been described in detail in previous reports, 
but basically a primitive form of L (J to) is first used which satisfies the tf>(w) 

bounds and then local corrections are made in order to push L(jw) towards the 
Tp(oj) , with iterations continued until the global 'error' (sum of distances from 
the i|/(w) ), is less than a prescribed amount. This is one of the most difficult 
parts of the program. 

A different approach was attempted, wherein equations have been developed 
which actually give L(joi) explicitly as functions of the bounds ift(w) . 

As expected, the equations are integral equations, so in the general case an 
iteration process is necessary to solve them. In the special case when the 
bounds are straight lines, no Iteration is necessary. Hence, one approach 

is to first use present techniques to obtain a reasonable first approximation 
to l-Qpt^^ > and then represent ^(w) as straight lines in the intervals where 
the approximate L o p t (jw) crosses the t// ( to) . As of this date, there has not 
been acquired sufficient experience to justify incorporation of either of these 
techniques into the automatic computer program. The analytic background and 
development is given in Appendix 1. 

2. Nonminimum-phase (nmp) systems . 

The computer program delivered Is confined to minimum-phase (mp) systems, 
i.e. for all possible plant parameter combinations over their range of uncertainty, 
the resulting plant transfer function P(s) has no right half-p’ane zeros (rhpz) . 
Then almost any set of performance specifications over almost any bounded parameter 
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uncertainty set (see (1) for limitations), are achievable. However, if P(s) 
has a rhpz, say at b then L(jw) = P(jto) G (jto) cannot achieve a crossover 
frequency (defined by f L (jto^) [ = 1 ) significantly more than b/2 (see 

(1)). Sensitivity reduction is achieved only for to < to^/2 approximately, so 
the adaptive capabilities of nmp systems are severely limited. Given an uncer- 
tain plant which is nmp for one or more parameter combinations and specific 
performance specifications, it is generally impossible (at least at the present 
state of the subject) to know apriori whether the problem is solvable, i.e. 
whether the performance specifications are achievable. One may therefore have 
to proceed with the synthesis details and be stopped (if the problem is not 
solvable), when he tries to find L t (j<o) — it does not exist. It is then 
necessary to relax the performance specifications in order to make the problem 
solvable. 

The objective was to Incorporate nmp systems into the automatic design 
program. However, it was first necessary to solve tne problem analytically and 
thoroughly understand the detailed steps In solving it, if done by hand or 
interactively with the computer — especially the reasoning involved in relaxing 
the performance specifications the minimum amount, should the attempted ones 
prove not solvable. Also, the proof that was used to derive the properties of 
the optimum L(ju>) in mp systems, was not valid for nmp systems, it was 
necessary to first establish on a rigorous basis the theoretical foundations 
of nmp systems. This proved to be a first class research problem in itself, 
as can be seen from the results achieved, described in Appendix 2. 


1 » 

3. Translation of time-domain bounds into equivalent frequency- domain bounds 

The synthesis technique operates in the frequency (w) domain, whereas perfor- 
mance bounds are often given in the time domain. The present program, therefore, 
has a routine for translating the latter into tu-domain bounds, but which does not 
work satisfactorily for higher order systems, where the effective 'time-delay' 

(say the time t ^ , for the unit step response to reach 0.1) is relatively large. 
The problem of translation from time to frequency is also more difficult in non- 
minimum-phase systems. Considerable effort was invested in obtaining a good under- 
standing of the time-frequency relations, in order to make significant improvements 
in this translation problem. The results are given In Appendix 3. 
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APPENDIX 1 

Explicit Solution for Optimum L ( j cu) 

There is given a boundary ip (u>) in the complex plane at each w , such 

that the loop transmission L(jto) must lie on or outside . A typical 

family of ip (to ) in the logarithmic complex plane (Nichols chart NC) where the 

abscissa B represents Arg L(jw) and A the ordinate represents f L ( jco) f 

in logarithmic units, is shown in Fig. 1. The reasons for the portions labelled 

, T in Fig. 1, have been given in detail in Reference 2. A certain excess 

e of poles over zeros is a priori assigned to L(s) , so that as s -*■ ® , 

L(s) ■+ — . The objective is to find L(jo)) which satisfies the boundary 
s 

constraints with a minimum value of k . Reference 2 details the reasons for 
this definition and the a priori assumption of a value for e . 

The ip(a)) are of two types. The one at w = ,5 in Fig, 1 is an example 
of the first which is single-valued for all B and exists for B from -360° 
to 0°, whereas ip^ i 5 an example of the second which is multiple-valued and 
generally exists for only part of B . It is assumed here that each ip has 
a continuous derivative at each point. Any sharp corners of ip can always be 
rounded with inf ini tessimal effect on ip , Consider a segment of ip in Fig. 2. 
At each point on ip a tangent line T^ can be drawn as shown at point P for 
example, whose equation is 

A cos 0 - B sin 9 = C (1) 

0 and will, in general, be functions of the point on ip , i.e. of the 

value of B being considered. The intercepts of this line are: on the B 

axis, B = -C /sin 0 and on the A axis, A = C /cos 0 . Suppose that the 
o o 
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acceptable In L(jto) at B » must be on or above B so that point Pj 
is acceptable. Draw a line Tj parallel to T Q , through P ] . The equation 
of this line is 


A cos 0 - B sin 0 = C 


( 2 ) 


with C > C . Let In L(jw) = A(to) + j B (to) , and suppose the boundary segment 
— o 

, . The condition that at B = B , A(w ) must 
1 P ' 


shown in Fig. 2 is for to = to 


lie above P , is expressed by writing 


A(B ,w.) cos 0 (B ,to. ) - B (to,) sin 0 ( B ,to ) = C > C (B ,to.) 
pi pl pi P J o p i 


(3) 


Note that C is a function both of B and of to as is 0 . The permissible 
o 

A is consequently also a function of these two, explaining (3). 

Consider the situation at point Q on i|>(Wj) in Fig. 2. We want 0 to 
be a continuous function of B , so clearly < 0 (B^) < 0 and C q (B^) > 0 , 
order that the A intercept is positive. The condition that L(jto ) lies 
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on or above Q if B = Im In L(jto^) = B^ , is the same as (3) if P is replaced 
by (1 , i.e. it is C ^C Q (B^,to^) . It is found that this is the condition at 
all points on i|>(Wj) ; also, it makes no difference if all or some of the tangent 
lines have A intercepts which are negative. Thus, it is noted that for $ 
of the first type, (i.e., 1 6 1 < ir/2^cos 0 > 0) , for all B values the accep- 

table condition is C > C 

— o 

Consider a boundary of type 2 with the acceptable region for L(jto) outside 
^ , as in Fig. 3* if we start at point Pj , define its 0^ > 0 and move clock- 
wise then the previous paragraph shows that from P| to P^ , the condition 
C >_ C Q applies. At ? 2 , 0 2 = -tt/ 2 ’ so Eq. (l) becomes B = C q . The acceptable 
condition is hence again C C q . At P^ , continuity of 0 gives, -tt < 0^ < 
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-tt/ 2 , cos 8^ < 0 and the intercept of the tangent line being C q /cos 0^ , the 
condition that L lie on or below P, again requires C > C because C/cos 0„ 

< C /cos 0 » i -F C > C and cos 0_ < 0 . At P. , continuity of 0 gives 

O j O j . n. 

-1 . Sir < 0^ < -ir , cos 0^ < 0 , so again C _> C q is needed so that L lies on 
or below . Thus, if 0 and C q are properly defined as functions of 6 

(and, of course, of to) the condition that L lies on or outsido a closed boun- 
dary is given by C >_ C q . 

The problem of choosing an optimum L(jto) can therefore be stated as follows. 
There is given an infinite set of equations, one for each to in [0,“) 

A cos 0 (to,B) - B sin 0 (to, B) = C Q (to,B) (4). 

A 

0, C q are continuous functions of to, B . Find In L(jto) = A- (to) + jB(to) , 
regular in the right half-plane and the boundary (rrhp), such that when the 
latter A,B are substituted into (4), they satisfy the equation with C >_ .C . 
Given that L(s) has a specified excess e of poles over zeros, so that 
L(s) -s- k^/s 8 as s + “ , the optimum L is that with minimum k^ . 

The requirement that In L(jto) - A(to) + jB(w) satisfy (4) with replaced 

by C •> C q , may be put in the form 

J^ e e j 9 (A + jB) = C (5) 

The general complex equation then has the form A + jB = e J (C + j D ) . One 
difficulty is that In e ^ Is not rrhp, because e ^ is all-pass. We 
therefore deliberately introduce e 0 '^ such that e a+ J® and j ts In are rrhp. 
This is done by choosing cr(to) as the Hilbert image of 0 (to) , using the Bode 
integral [1] 










A] -4 


2 /CO 

- a(o) = 2£L 

T J 0 


9 (y)dy 

, 2 2s 

y(to - y ) 


The desired equation is now 


e^' 6 (A + j B) = (C + jD) e° = C + jD 


whose real part is Eq. (5) multiplied by e In order to ensure that A + jB 

A -S a 

is rrhp as required, D must be the Hilbert image of C = C e , i.e. using the 
Bode integral [1], 


D(w) e 


-2(0 f C (y) e a ^ dy 
^ J 0 to 2 - y 2 


Then (7) can be written as: 


A((o) = C(go) cos 0 (w) 


2hie~ 0 ^^s?n 9(h)) f“ e 0 ^C{y)dy 


n 2 2 

0 (i) - y 


BM - -C(B) sin eW -2j5^2LlM f“ e° ( ^C(y)dy 

n J 0 (o 2 - y 


It is understood that in the general case A, C, 0, etc. are functions of B 

as well. Recall that the objective is to minimize k_ where 1 im L(s) = k 0 /s e . 

** s-* 33 * 

Hence, lim 1 n [ L (jtu) | = lim A(w) = In L - e In w and In k- = lim[e Inw+ A(io)] 
(0 -x° oj-^ 30 

Minimization of k^ means minimization of lim A(w) . In Fig. 1, by recalling 

the evaluation of 0 by means of the 'continuity of G , it is seen that 

lim 0(w) = -ir/2 , so the first term on the right side of (9a) disappears and the 
(0-*=° 


%3f. |f i ; -m. 


I ; M"". : -ytr 


" Vc • -f.. 


AT -5 


second term becomes 


2ue 


-ct(oj) 


-eo 


a CT ( y ) C(y) dy 


J 0 


2 2 
0) - y 


( 10 ) 


whose limit as to « is to be minimized. Since every factor in (10) is 
positive (recall cr(y) is real) except possibly for the free C(y) , it is 
clear that a necessary condition for minimization is to make C(y) as small 
as possible. But C(y) C Q (y) . Hence, choose C (y) = C Q (y) . This proves 

that the optimum L(jm) (if it exists) lies on its boundary (to) at each tu . 
It has previously been proven [2] that if any solution exists at all to the 
problem, then an optimum must also exist. If the plant transfer function set 
is minimum-phase over its entire range of uncertainty then it is known [3] that 
a solution exists. 


Application to Numerical Work 


The simplest application of the above under consideration for numerical work, 
is to first obtain an approximation to L q by means of the computer program 
previously delivered to NASA. This approximation gives one a good idea of the 
region on each 4>(w) , where L ^ crosses. A straight line approximation, is 
then made of this region. Then 0 (w) and C (w) = C o (w) ’ n (9 a »b) are known 
as functions of w , and (6) can be used to evaluate a (to) - a(0) . Note in 
(9a, b) that the constant o(0) cancels out and need not be known, so it can 
be set at zero. On the vertical line of the 'universal high-frequency boundary 1 
in Fig. 1, C = B a constant, corresponding to the value of the constant 
phase. It is clear from Eqs. (9a, b) that the straight line approximation is 
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the simplest to use because no iteration Is necessary as C q , 0 are then not 
functions of B . Thus, A .and B can be calculated independently. 

Experimentation with this approach has not as yet been completed. In 
numerical examples where the are well approximated over a significant 

interval near the point where L A p PROX (j“) intersects (f; (to) , the results have 
been quite good, permitting quick convergence to the optimum L . However, it 
must be noted that the final result gives In L(jtd) = A (to) + jB(to) in the form 
of numerical data, not as a rational function. It is necessary to use the RFA 
(rational function approximation) subroutine ot the computer program package 
previously delivered to NASA, in order to secure L(s) as a rational function 
approximation of the numerical data. In numerical examples where the straight 
line approximation is valid only over a narrow interval in the region where 
*'APPROX^ 0l O intersects i|j(w) , the results have been poorer, as could be expected 
Application of (9a, b) still gives a clear-cut answer, but on some <J/(co) the 
results put the point on a part of the straight line where it is no longer a 
reasonably accurate representation of ^ (to) . One must then take a compromise 
straight line in between, i.e. iterate. If so, one might as well incorporate 
the iteration directly Into (9a, b) by using a second-order approximation of 
iK<i») in the region of Interest, at least at those to values where it appears 
necessary. Then 0 and C — in (9a, b) must be written as functions of B 
as well as of to . There has not been sufficient time, as yet, to thoroughly 
evaluate these alternatives and incorporate them into the computer program. 
However, we are confident that this approach will significantly improve the 
program for finding the optimum L , and intend to continue working on this 
phase. The improved program will be made available to NASA when it is ready. 
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APPENDIX 2 


OPTIMUM SYNTHESIS OF NONH I N 1 MUM-PHASE 
FEEDBACK SYSTEMS WITH PLANT UNCERTAINTY 


Isaac Horowitz* and Marcel Sidi* 


ABSTRACT 


Linear, time invariant feedback systems in which the constrained plant 
transfer function has right half-plane zeros, are perforce nonminimum-phase, 
and their attainable benefits of feedback are inherently restricted. This 
paper presents criteria for determining whether a given set of performance 
specifications are achievable, and if so, a synthesis procedure is included 
for deriving the optimum design, defined as that with effectively minimum loop 
transmission bandwidth. The properties of the optimum design are derived and 
its uniqueness proven, for both the minimum and nonminimum-phase feedback 
systems . 


*Department of Applied Mathematics, The Weizmann Institute of Science, Rehovot, 
Israel, and the Department of Electrical Engineering, University of Colorado, 
8oulder. 

This research was supported in part by the U.S. Air Force Office of Scientific 
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Administration under Grant NSG 11^0, at the University of Colorado. 


OPTIMUM SYNTHESIS OF NONMINIMUM-PHASE FEEDBACK SYSTEMS WITH 
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PARAMETER UNCERTAINTY 


I. INTRODUCTION 


There is given a linear time- invariant constrained plant 
with minimum-phase (mp) transfer function P(s) (see Definitions II.). 
There is uncertainty in the parameters of P(s) , but the ranges of 
the uncertain parameters are known. The uncertainty may include the 
excess (e^) poles over zeros of the mp plant, up to three Cl], 

e.g. H 6^ n , and parameter ranges in which the plant is 
unstable [ 2 ]. It is desire,d that in response to command inputs 
r(t) and disturbance inputs d(t) in Figure 1 , the system output 
should satisfy, over the entire range of parameter uncertainty, 
specifications of the form 


A t (u>) S |T( jw) I £ B t (w) 


T ( s ) 


A c(s) 

" RTsT 


A^Cw) « |T d (jui)| £ B^ (to ) 


T d (s) 


£ C( s ) 
D(s) 


( la,b) 


for given A^tto) S B d (to) . 

It is possible to demand A^(w) = 6^(10) at a finite number of 
( 1 ) values, at which the loop transmission must then be infinite. 
Although these performance tolerances are in the frequency domain, 
it has been shown [ 3 ] that they suffice to satisfy time-domain 
specifications of the form a^ m (t) * c^ m ^(t) £ b^ m (t) i = r,d 


■' /) 2_2 

where c| ra ^(t) is the mth derivative of c. (t) for m = n 

finite 3 and a^(t) = fc^(t) allowed, at most, only at zero and 

infinite t . 

* 

A detailed synthesis procedure has been presented [2] for the 
above problem, and an existence theorem [13 guaranteeing its 
validity, providing PCs) is mp over its entire range of uncertainty. 
In addition, under a certain definition of optimality and under 
simplifying restrictions reasonable only for a mp plant with a 
relatively large range of uncertainty and/or narrow performance 
tolerances, the properties of the optimum design were derived and 
its existence and uniqueness proven [4]. This paper extends the 
synthesis procedure and the optimization results to both mp and 
nonminimum-phase (nmp) plants in which most of the simplifying 
constraints are eliminated. It is necessary in any case to remove 
the simplifying restrictions in order to cope with the nmp system. 
Additional reasons for removing them are: 

1. They do not necessarily hold in problems with moderate parameter' 
uncertainty and/or large performance tolerances. 

2. It has been shown Cl, 3] that the above linear, time-invar5i\H t 
(&ti) synthesis technique based on frequency response, may be 
rigorously used in linear and nonlinear time-varying systems with 
uncertainty. The key ingredient is the replacement of the non- 
linear (or linear) time-varying plant set by an equivalent £ti plant 
set. The restrictions which reasonably apply to a true £ti plant 
set (i.e. from a S-ti model) may not at all apply to this equivalent 
&ti plant set. Also, this equivalent £ti plant set, obtained from 

a nonlinear or linear time-varying plant set, is much more likely 
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to be nmp. This is one of the principal motivations for solving 
the nmp problem - to greatly extend the £,tv and nonlinear problems 
solvable by the ’equivalent £ti' method. 

3. Optimum &ti multiloop and multivariable synthesis based on 
frequency response, is to a large extent a matter of the optimum 
shaping of analytic functions (the loop transmissions) in the 
complex plane. While some progress has been made [8,9], the 
optimization problems are in general unsolved. The techniques and 
theorems used in this paper for extending the single-loop results, 
are useful and suggestive for the more difficult multiloop problems. 

It is noted that the nonminimum~pha.se concept -ms first 
introduced by Bode [7]. The problem in control systems has been 
treated for deterministic and stochastic inputs, notably by 
Chang [14], Horowitz [5] and very thoroughly by Youla et al [15], 
However, these treatments have been confined to systems with 
no parameter uncertainty. 
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I i. Definitions . 

Throughout this paper, the left half-plane (s,hp) includes the 

imaginary axis, while the right half-plane (rhp) excludes it. If 

F(s) is a rational function in the complex variable s then it 

is here called minimum-phase (mp) if it has no zeros in the finite 

rhp. Otherwise, it is nonminimum-phase Cnmp). A simple way to 

handle a nonrational function is to approximate it by a rational 

function over the frequency range of Interest and deal with the 

approximation. In this way it is seen that any reasonable 

-sT 

approximation of e is nmp. Note that rhp poles are admissible 

in a mp function. 

This paper is concerned only with problems in which the 
relevant functions can be approximated as closely as desired by 
rational functions, so the above definition suffices. 


* * * 
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I 2. Constraints on LCs) , T(s) due to 


n _ m P_ 


PCs) 


Figure 1 presents a two-degree-of -f reedom [5] structure 

in which F(s), GCs) are to be chosen to satisfy the system 

'ji-'U-'L.- 

specifications Cl). If PCs) has one or more rhp h at 
let them be explicitly shown by writing 




q 


• PCs) = P, Cs) Cl-T-. s) ... (1-T s) 

11 q 

= <P 1 (s) 3 (l+T iS )>[ ={P m Cs)>[A(s)] , (2a, b) 

LCs) = GCs) P(s) = (GP )A(s) * L Cs) ACs) . 

m m 

PCs) , L Cs) are mp while ACs) is an 'all-pass' function, 
m m r 

because }A(jw)| = 1 for all w . If the system is to be stable 

then it is impossible to cancel the rhp zeros of P by means 

of rhp poles of G because cancellation, being imperfect, results 

in rhp dipoles for the loop transmission L = GP and therefore 

for T = FL/Cl+L) . Hence, both L and T are constrained to 

have the rhp zeros of P as zeros. While this is actually also 

true for the £hp zeros of P , these can at least be cancelled 

at a nominal set of parameter values , so the resulting nominal L Q , 

T q can be free of these zeros. In the nmp case they must appear 

even in the nominal L , T . 

o o 

Consider the effect of a single^ all-pass section 

ACs) = Cl-xs) / CI+ ts ) in L(s) = L (s) ACs) . The phase of ACjw) 

m 

O 

decreasess monotonically from zero at m = 0 , to -90 at m = 1 /t 

o 

and tends to -180 as co ■* m . ’ In a realistic feedback design. 


fit-*. 

|L(jw)j >1 by a reasonable margin over some w range, in 
order to obtain the desired benefits of feedback. Beyond this 
range it is important to reduce |L(jm)| vs w , as rapidly as 
possible [ 2,8]. One pays with phase lag for the rate of reduction 
of |L(jtu)| , with the stability requirement limiting the fund of 
phase lag available. ACjw) contributes phase lag, thereby 
diminishing the fund, without contributing anything, being all- 
pass, to the reduction of |L(jw)| . At w = 0.5/r, the 

o o 

amount of phase-lag lost is 53.2 which when added to say 35 

o 

for phase margin, leaves only about 90 for magnitude reduction. 

It is therefore impossible to obtain a crossover frequency u> c 
(defined by |L(jaj ) | 1) , much larger than 0.5/t . The above 

discussion is qualitative. The quantitative limitations will be 
forthcoming from the quantitative synthesis procedure in this paper. 

XI. SYNTHESIS PROCEDURE FOR THE NMP SYSTEM 

The synthesis procedure parallels that for the mp system 
[2], with certain modifications. The frequency response 
specifications are assumed to be given by (la,b). The problem 
of translating time-domain into w -dom .in tolerances in mp 
system has been discussed [2,6]. While no rigorous translation 
exists, there is in practice little difficulty in effecting one 
satisfactory for engineering purposes. Appendix 5 discusses 
the additional complexities that arise in the translation, due to 
the system being nmp. In presenting the technique, primary 


6 . 

attention will be paid to (la) (see Figures 3, 4 for a typical 
example of (la)). It will then be obvious what to do about (lb). 
Also, to simplify the presentation, the order of the plant is 
assumed known. The modification needed for a plant of uncertain 
order (up to a factor of three) will then be evident. 


Let all the rhp zeros of the plant and only these, be 
contained in the polynomial N(-s), -s being used to emphasize 
its nmp character. Let N (-s) be the nominal value of N(-s) 
obtained at a set of plant parameter values arbitrarily defined 
as nominal. There is therefore no uncertainty in N Q (“S). P(s) 
may also contain poles whose range of uncertainty includes 
portions of both the &hp and rhp . Denote these by the 
polynomial d(s), with nominal value d^(s) deliberately chosen 
with all its zeros in the Jihp. 


The above 
by writing P 

L = PG = 


fhp zeros and poles are explicitly displayed 
= N(-s) P-^(s)/d(s) , and 


N(-s)P 1 G 

■N 0 C- S) - 


[■p in (s)N (s)Gl 
10 o 


[N(-s)P,(s) d (s) I 
1 o 

d(s ) 

U Cs) J 


d (s) 

A 


- N ( -s ) P, n (s) d(s) ■* 


= CA o (s)][L mo (s)]EV(s)] = L no V(s) 

where P^ Q (s) is the nominal (and hence not uncertain) represent- 
ation of P 1 (s) when the plant parameter have their nominal values. 
V(s) is the only unceratin part of L and is unity at the 
nominal values of the plant parameters. L mo (s) is mp and stable 
of course, and L nQ is the nominal nmp loop transmission. Since 
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T = F nc > 

A In L = i Jin V 
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(4a,b,c) 


due to uncertainty. As the plant parameters range over their 

regions of uncertainty, A Jin L(jw) = A Jin V(jw) occupies a 

known region in the complex plane at any fixed co . It is a 

straightforward matter to use (4b) to find those values of 

L (jw) = A (j<jj) L (jw) which satisfy (la) over the entire 
no o mo J J 

range of A Jin V. The boundary of the acceptable L nQ is 

n i 

denoted by 3 , and of the corresponding L - A - L by B 

to mo o no ^ to 

The logarithmic complex plane (Nichols chart) is useful for hand 
calculations and for seeing the design 'tradeoffs' [2, 8], but 
of course the procedure easily lends itself to computer calcualtions . 
The procedure is illustrated by the following simple example. 


Example ♦ 


PCs) 


k(l-TS ) 

s(l+|3s) 


with parameter uncertainty: 


k e [i,3] , 


B £ [0.3,1] , t 6 [.05,.l] . The arbitrarily chosen nominal 
values are k Q = 1 , B q = 0.3 , t q = .05 giving 


- k( 1-ts ) ( 1+ . 3s ) 

V(s) = 

(1-. 05s) (1+Bs) 

The range of V(j4) on the Nichols chart shown in Figure 2. Note 
that it includes the point 1 = 0 db u . Suppose that }T(j4)j 
is allowed to be between 0.5 and -3.5 db i.e. A Jtn | L€ j 4 ) J £ 4db 


Using (4b) it is quite easy to find the permissible values of 


L mo (j 4 ) of (3), such that A An|T(j4)| S 4db for all V(j4) , 
a portion of which labelled as is shown in Figure 2, i.e. 

A £n|T(j4)| 4db so long as L nQ (j4) lies on or above . 

Such bounds on the nmp L nQ exist for each to £ [0,«) . How- 
ever, in actual numerical design it is more convenient to work 
with a mp function (the Bode integrals [5, 7] can then be used), 
so the corresponding boundary of L^Cjw) of (3) is easily 
derived as that of L nQ (j4) multiplied by A~ i (j4) = /22 . 6 , 

in this case, giving the boundary B^ in Figure 2. Similar 
boundaries exist at each to , with the amount of shift in angle 
due to A (ju) , tending to 180 as to becomes large. 

IX 1 Disturbances and the Single H ' sh-Frequency Boundary 

In the above, only (la) was considered. The disturbance 
specifications (lb) may require all or part of the boundary to 
be changed, so as to satisfy both (la,b) . It is important to note 
however, that even in the absence of significant disturbances 
(so assumed in the balance of this paper), the latter should 
not be entirely ignored. If it is, then design economy (as will 
be seen) will result in very large peak values of |L/(1+L)| at 
higher to . Such peaking does not violate (la) because in (4a), 
the prefilter F is abailable to drastically attenuate the 
higher frequency components of the command input. It is however, 
assumed impossible to prefilter the disturbance D in Figure 1. 
Hence, one should always add a specification of the form 


\ 
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s Y 


(5) 


with y chosen by the designer. This leads eventually, in the 
S,ti mp problem, to a single high frequency boundary , 

effective for all to greater than some to^ , for the following 
reason. 

e 

At large to , P -+■ k/s ^ so VCs) -*• k/k Q whose range 

( -template .) is a vertical line in the Nichols chart, of magnitude 

20 logCk /k . )db . Assuming y in (5) is a constant in this 
max mxn 1 

range, the resulting boundary of L (jw) is the sarae f° r 3JL1 

to for which the above is applicable, e.g. in Figure 2 

for y = 2 . 3 db , k /k . = 10 . In nmp systems this 'universal' 

1 ’ max mxn c J 

high-frequency boundary of L nQ is translated to the right by 

e 

|Arg A q C j ui ) j in ordr to obtain the boundary of the mp L mQ (j(o) . 

The resulting boundaries are denoted by B. to emphasize their 

dependence on w . Figure 2 displays some for 

A (s) = (1- . 05s ) / ( 1+ . 05s ) . Note how the corner point V 
o w 

moves to the right as w increases. In general, Arg A o (jm) 

O O 

tends at laest to 180 (if A (s) is of order 1, 180n for 

o 

order n) , so that sooner or later the right side of B h[D lies 

in the positive angle region. It is this excursion of B hu) 

into the positive region which may render it impossible to 

satisfy a given set of specifications. It is proven in Appendix 

4 that if L (iw) cannot reach the corner V before V 
mo J to tu 

is on the zero angle line (i.e. while Arg L is still 
negative) , then it is impossible for L mo to satisfy the 
specifications and also 0 as tu ro . 


f\z~ 10 - 
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II 2 Design Optimization 

In mp systems, any realistic specifications of the (la,b) 

type (see Cl] for precise conditions) , may be satisfied for a 

large class of plant parameter uncertainty problems, by an 

infinitude of designs. A nmp problem if it is solvable at all, 

will also generally have an infinitude of acceptable designs. 

Uniqueness is achieved by defining an optimal design, here chosen 

as follows. Let e be an apriori specified excess of poles over 

zeros for the nominal loop transmission L q (s) j so that as 

0 

tj co , L (s) -j-k/s . 

5 o o 

Definition : The optimal design is defined as that which satisfies 

the system specifications with minimum value of k Q . 

This definition is motivated by the practical importance of 

decreasing |L(ju)?| vs w as fast as possible, in order to reduce 

the effect of high-frequency plant parasites and of sensor noise 

[2, 5, 8]. It is noted that essentially this criterion was used 

by Bode [5, 7] for the derivation of his 'approximately ideal 

characteristic.' It is important to note that by working with the 

mp L (s) of (3) instead of the nmp L (s) , the problem of 

finding L Qpt becomes one of finding (L mo ) t . Henceforth, 

L ^ denotes (L ) . , unless otherwise specified, 

opt mo opt 

The following property of L Qpt is worth noting, derived 
from the Bode integral [5,7] 

00 

&n| L( 0) | - Jln| L( ro )| = - - / Arg LCjw)d£n u . (6) 

— oo 

This relation may be applied to our problem by modifying L at 


Dz- u. 


very large w , to be nonzero as w -> °° and, if necessary, at 
very small w to be noninfinite as w ->• 0 . It shows that 
the optimum design strives to maximize its phase lag, subject 
to the boundaries B 

to 

It is proven in Appendices 1-3 that if any design exists 
which satisfies the specifications, then L 0 p-j- is unique and 
exists in a limiting sense, i.e. may be approximated as clearly 
, as desired. Also, L ^ lies on its boundary ar every to 

from zero to infinity. It is crucial however, to carefully 
include certain conditions in order that the problem be well 
defined. Figure 2 may be used for this purpose. Let 
be the applicable 'universal' high-frequency boundary in a 
purely mp problem, i.e. A q (s) = 1 . As to ->• «* , L 0 pp must 
at least asymptotically approach the vertical line W" YZ 

O 

(in Figure 2) corresponding to Arg L - -90 e . L ^ must 

therefore make the transition from B, to this vertical line and 

h 

as so:n as it 'turns the corner’ at QV there is nothing to 
stop it from very rapidly increasing its phase lag. However, a 
rapid increase in phase lag is accompanied by a sharp increase 
in magnitude [5, 7] so L would follow QVW and then 

instantaneously increase its phase lag to -90 e degrees, where- 
in | k Q p.i_ | abruptly goes to « along QVWW' and returns on 
W"YZ . To prevent such a discontunuity , and for obvious 
practical reasons, the supplementary ’transition 1 boundary W-jjV^YZ 
is added [4], i.e. that for all w more than some , L must 

lie in the closed region bounded by ..QVW W. V^YZ . 


The above suffices for mp Ati problems with large parameter 

uncertainty and/or narrow performance tolerances. In these, 

L opt^“) must be so large over a large to range, that it can 

reach the NQV region in Figure 2 only at w values so large 

that the plant function P(jw) is very close to its asymptotic 

-e 

high frequency value ks p . The ran^e (template) of V(jw) 
of (3, 4) is then a line so is the single effective universal 

boundary for all higher w . This is not necessarily so for the 
problems listed in the Introduction. In these L nQ may cross the 

o 

-180 line before P(jw) has degenerated into ks F . The 
corresponding B^’ need not form a nested set in this frequency 
range, but instead cross each other, expand and contract in size 
and shape etc. A more careful definition must be made of the 
transition from the boundaries B n obtained from Equation (4). 

O 

to the high-frequency asymptote Arg L no (jw) = -90e , in order 

that a continuous L 0 pt^ w ^ may resu l t * To do this note that there 

must exist aw , 3 , V w > w , the boundaries of L do 

3 no 

not enclose the origin of the arithmetic complex plane, and lie 

o o 

entirely in the second and third quadrants, -270 < Arg < -180 

Otherwise, it is impossible for L^Cjw) to approach zero as, 

O 

m -»• oo . Sooner or later L (qw) must cross the -180 line 

no 

between the origin and the bouriaries , and some interval after 
it has done so, its phase may be sharply decreased. seeks 

to ’escape' from the boundaries to its final alymptote as quickly 
as possible. The transition boundary consists of three parts. 

o 

The last is a segment of the line Arg L no = -90e (e.g. ODE 
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in Figure A5a for the case e = 4) , corresponding to YZ of 
Figure 2. The second is one or more portions of an arc of constant 
magnitude |L no | = CM corresponding to W^W-^ of Figure 2 
(see Figure A4 for CM > 1 and Figure A5 for CM < 1) . 

THe first part (p^) of the transition boundary must join 
CM to the B n and is defined as follows. For each B n , w > to" , 
choose a straight line tangent to B™ and terminating on CM 

(e.g. in Figure A4a where B^ = B^) . If several such points 

of tangency exist, choose one for which p does not intersect the 
interior of B^ . If CM intersects B™ , p y may, of course, 
not exist. Henceforth refers to the boundary modified by 

and it is therefore also necessary that the p be chosen so that 

0J 

this set {B n } of modified boundaries continue to satisfy (Al)- 
It is conceivable for L no (jm) to be on CM for some to range, 
be forced to return to B n , but after one or more such excursions 

U) 

eventually return to CM while in the second quadrant, and leave it 

o 

then only for Arg L =«>9 0e 

® no 

Given the above problem definition, it is proven in /ppendix 

1 that L . must lie on its associated B n at every to , in 
opt to 3 * 

Appendix 2 . that if any L exists which satisfies the B^ then 
L opt exists 
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exists, it is unique. The 


and in Appendix 3 that if L ^ 
above apply to both mp and nmp systems. It is well known and 
formally proven in [2] that in the mp case a solution does 
exist. It is possible, however, that no solution exists for 
a specific nmp problem. Appendix 4 presents some criteria 
for determining whether this is so. 


III. DESIGN EXAMPLES 


Example 1 

The plant transfer function 

p (s) = with k e [1,3] , b £ [0.3,1] d £ [0.05,0.1], 

The step response bounds were originally those shown as I in 
Figure 3a, resulting in the equivalent frequency-response 
bounds on the magnitude of the minimum-phase T^Cjw) (see 
Appendix 5) I in Figure 3b. The disturbance response parameter 
Y of Equation (5) is 6 db . 

Design 

The range (templates) of V of Equation (3) were obtained 
and used to find first, the bounds on L and then as 
discussed in Section II (Figure 2) , the resulting bounds B u 
on ^ (jw) > shown in Figure 4a. Using the technique of 
Appendix 4, it is found that L mQ cannot be reduced sufficiently 
by the time Arg L is zero degrees. Therefore the specifications 
on the system step response must be relaxed. In Figure 4a 


p\Z- 15. 

it is seen that primarily the bounds for to € [2,10] must be 
relaxed as B 2 compels |L mQ (j2)| to be^ 0 db and compels 

Arg L mo ( j 1 0 ) to be rather large. A new set (II) of frequency 
domain bounds and their corresponding time domain bounds are 
shown in Figures 3b, a respectively. The resulting new bounds 
on k mo £jio) are shown in Figure 4b together with the L mQ (j(o) 
used. Experimental design results are shown in Figure 3a. 

Example 2 . 

In this example P = ( 1-ts )ka/s( s+a) with uncertainties 

k £ [1,10] , a £ [1,10] , and two separate values of t were 

considered: 1) t = 1/4 2). T = 1/60 . The bounds on the nmp 

’ shown in Figures 5a, b^ are therefore the same in both 

cases. But the bounds on the mp for T = 1/4 in Figure 

5aj are shifted to the right much more than those in Figure 5b, where 

r = 1/60 . A computer program was prepared which derives an L(jw) 

function which lies on the boundaries [9]. In Figure 5a this 

program perforce gives an L function which crosses into the 

positive phase region, as shown. Such a L function cannot go 

to zero at large u . In Figure 5b, L is attainable, and the 

mo 

resulting nmp L no (jm) is also sketched. With such a computer 
program it is a simple matter to check whether any specific nmp 
problem is solvable. 


I 


IV. CONCLUSIONS 


Optimum single-loop Jlti synthesis for plant? with parameter 
uncertainty, previously developed for a special class of mp 
problems [4], has been extended to a much larger class of mp 
systems and to nmp systems as well. This extension permits 
the 'equivalent Jlti' technique developed for linear time- 
varying uncertain systems which are effectively mp [1], to be 
also applied to JLtv uncertain systems which are effectively nmp , 
the resutls to be made available in the near furture. The results 
of this paper are also obviously applicable to Jlti synthesis 
techniques previously developed for more complex mp structure 
such as the cascaded multiple-loop [8] and the multivariable [9]. 
In all of these, the availability of a computer program for 
finding the L(jto) function which satisfies the bounds, makes 
it an easy matter to determine whether the original nmp 
specifications can be satisfied or whether they must be modified. 
The synthesis technique provides insight into the minimum required 
modification. 
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APPENDIX 1. 

L opt (ja>) LIES ON ITS BOUNDARY 

It was shown in Section II that the problem of optimization 

of the nmp L nQ (s) of (3), may be transformed into that of the 

mp L (s) . Until specified otherwise L , L . , L , B etc 
mo ? opt a w 

refer to the mp items. The following restrictions on the 
boundaries are highly realistic and obvious in any real 

problem, and are listed for the sake of a clear mathematical' 
problem statement: 

bound simply connected regions in the complex 
plane. They are continuous on [0,«>) in the sense 
that given any point z on , then for Aw small 

enough, there is a point z+Az on B^^ is 

as close as desired to z , in the Euclidean metric. 

Also each B is a continuously differentiable curve 
in the complex plane. 

The continuity requirement does not require similar continuity in 

(D of the range (template) of P(jw) . Even if the original 

is not continuously differentiable it can easily be deformed to 

be so with infinitesimal changes in its actual shape. The region 

o 

bounded by B which is forbidden to L(iw) is denoted by B 

J W ■ J (il 

Theorem 1. L„^(iw)if it exists lies on B 

opt ■’ U) 


A2 -)2 


Proof. Let L_(jcu) be submitted as candidate for optimum and 

- g 

suppose there is at least one to interval (not a point because 
of (Al) ), denoted by I = (A,B), in which L a j. itself. 

Suppose that it is always possible to v i , ind an H(jto) such that 
L H(jw) & B° for any to , with Jtn | H ( 00 ) | < 0 , and 1+L H 

a. /CO cl 

having no rhp zeros . Note that both L and L_H , as 
candidate functions, must have the assigned excess e of poles 
over zeros, so H(«>) must be finite. Then L a H is better than 
L a in the sense of the definition of optimality. The theorem 
is therefore proven if such an H(jfu) can be found whenever such 
an interval I = (A,B) exists. 

The main problem in choosing H is to ensure that 
L a H(jw) £ B° , V to. The problem does not exist for to £ I = (A,B) 
or in any other intervals in which L £ B , because L(jw) 
may thelr*e be shifted in any direction by a finite amount without 
entering B°(to) . Thus, the worst case to consider is when 
L (jto) £ B , to £ T (complement of I). Consider in Figure Al , 

a. (0 

a very general form for in the arithmetic complex plane, 

(o^ £ X . In Figure Al , suppose L a (j(o-^) has the value A £ B y 
If Arg H(j<o^) = 0 and Jin |H(j(o^)j £ 0 , then clearly 
L H(jw-,) £ B° . If L C jtu , ) = E in Figure Al , then 

ct X tO^ <3. X 

Arg K(jto^) = 0 , Jin | H(jto^) j 5 0 is satisfactory. This strategy 
(Arg H s 0 , Jin J H | positive or negative) fails only at points 
such as C , i.e. if L (joj-i) = C . Such points are temporarily 
ignored as if nonexistant. Let 

Jin H(jto) = Jin | H | + j Arg H = H^(w) + j H 2 (u>) 


( A2 ) 
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and choose Arg H^(to) = 0 , V w £ I . It is required that 

£ 0 in those intervals in I in which L a (jw) has the property 
of A in Figure Al, and that H 1 i 0 in those intervals in which 
L a (jw) has the property of D or E . These requirements dictate 
say, n zero crossings of H-^Cco) at in T = ( 0 ,A) U (B,«0 

(see Figure A2), giving the n+2 intervals = ( 0 ,to^) , . . . = 

(rn^jA) , ^+2 = ^ ,w k+l ^ * ,R n +2 = ^n’ 00 ^ • Only the values of 

HnCu) in I are available to secure these required properties for 
H-^Cm) , inasmuch as ^(w) has been chosen zero for w £ I . The 
candidate can be such that a combination of the following may 

be necessary for H : 

(Aa) H^(m) 5 0 in R^,Rg,... in (0,A), with H-^ (to ) € 0 in 
Rj ,R^ , . . . in ( 0 ,A) 

(Ab) the opposite of (Aa) , i.e. H^w) s 0 for R^,Rg 5 etc * 

(this case is shown in Figure A 2 ). 

(Ba) If the number of 'intervals’ in (B,®) is even then the 

firs, third, etc. of them must be > 0 and the others z 0 , 
in order that H^(®) < 0 . 

(Bb) If the numbers of 'intervals’ in (B,®) is odd then the 
first, third, ... must be £0 etc. 

The following combinations are possible: (1) Aa , Ba (2) Aa, 

Bb (3) Ab, Ba (4) Ab , Bb, and the objective is to prove that 

can be constructed to give any desired one of these four combinations. 

Obviously H-^(Q) = 0 may be assigned for all cases. 
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Consier the Hilber transform relating the real ( H 1 ) and imaginary 
(H 2 ) parts of the analytic function Jin H(s) , s = jm[5,p. 309], 


E* (w ) 


& H^ftO-H^O) 


2 

0) 


” H 2 (C)dC 

J o n 

0 c(u -S ) 



eu)dc 
2 ,2 

w -C 


(A3) 


with H 2 (£) = 0 deliberately picked V hi £ I , giving the limits 
A, B instead of (0,°°) . Equation (A3) is used to prove that 
the desired properties of in I are achievable by a suitable 

choice of H 2 in I . The technique is borrowed' from those used 
in Total Positivity [10, 11]. Pick a continuous H 2 (?) in I 
with n zero crossings A £ < n 2 ... < n n £ B and in addition 

H 2 (A) = H 2 (B) = 0 . 

/ W 1 » ,(ii n\ 

Step 1. Define K( _ as the n x n determinant whose 

— «i> v 


2 2 -I 

ij th element is k(w^,Cj) = (m^-i;.) 


H- 1 « S n- and 

J ^ J J 


to i £ I ordered as in Figure A2. It is now proven that K(...) has 
constant sign for any ok pair satisfying the above ordering 

conditions i.e. k(w,c) is totally positive [10]. Since each 
C i , oij can be continuously varied (so long as the ordering is 
kept), a sign change of K(...) requires K(...) = 0 at some 


, Wj pair. 


Let the i th column of K(...) be the vector 


If K(...) = 0 at some e. , w . pair, then 3 scalars a.,... a 

„ 1 J n in 

2 . 2,-1 


n 

not all zero, 3 I a. a. = 0 . Hence 

1 


I X - 


= 0 


for j = 1,2,..., n, which means that E = N(u) 2 >/ir(w^-c?) 


This is impossible 


has n zeros at u>. , j = l,2,...,n 

2 2 2 2 

because N(co ) is of degree n-1 in to and (to £ 0 


\ 
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for w € I , £ I , The only other possibility that K = 0 

is impossible, as follows. 

The determinant KC..) is expanded by its first column, 
n T 

U i 2 A 9 

giving K(..) = J — 2 2 = ^^l^ = t M 2 - L ) letting z 1 = with 

1 w i~^l 

J. the cofactor of the (i,l) element of the corresponding 

matrix, and r, 1 does not appear in any J. . Thus K(..) can be 

2 

regarded as a function of £, = z^ . If KC..) = 0 then the 
derivatives ifi^ m ^(z,) = 0 , m = 1,2,.. . if/^Cz,) can be 

1 . rf 1 

obtained by differentiation the first column of KC..) with 

respect to z^ , giving a determinant whose first column is 
2 2 

l/Cw^-z^) ; i = 1,2,.. n , while all other columns are the 

( JTl) 

same as in K( . . ) . Similarly if; 'kz^) consists of KC . . ) with 

2 m+1 

the first column only replaced by m! /Cuu-z^ . As m 

becomes sufficiently large, one term will dominate in the 

expansion of ij> m (z^) by the first column - the term with the 

smallest Iwf-z^l • Hence only its cofactor need be considered. 

This process is repeated, with the cofactor considered as a 
2 

function of > etc until the final cofactor to be considered 

2 2 — 

has '■nly one term l/(w.-C ) t 0 because u- E I , t E I . 

3 n j n 

Hence K( . . ) has the same sign V w- £ I , <;• E I , ordered 

^ J 

as previously indicated. 

Step 2. From A3 and the definition of kCoj.,?.) in KC-.) , 

B 3 

h1(w) = / kCm ,£ )Q ) dc , in which from (A4) below, 6(<;) has 

x A 

n+2 zeros 9t n = A < n. < ri„ . .. < B = rj , t i and has 

'o '12 n 'n+1 
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the same sign in between any two of these consecutive zeros 
and the opposite sign between the next two. It is next shown 

ft 

that of (A3) cannot have more than n internal zeros 

in I . Assume this is not so and that it has n+1 internal 
zeros at 0 < < w 2 . . . < tJ n+ ^ in T . Let io n+2 be any 

internal point in T , and consider the n+2 square matrix M 
whose ij th element is 


I 

"j-l 


kCtcu,?) |0(5)|dC 


for i = 1,2,..., n+2 ; j =1,2,... n+1 , 


while the (n+2) column consist of of (A2), i = l,...,n + 2 

(Det M ) is zero because from (A2,3), its last column is a 

► A 

4 + - 

linear combination of the others. Since by assumption, = 0 

for i = 1,2,..., n+1 , expansion of det M by its last column, 
gives 


n l B (u- 
I det M.j = | H"C« +2 )| f ... f K [ 1 
1 n 4 A r\ 4 1 


,w ii 

’ n+1 




n+1 


n+1 

, n I d ^i* * * d? n +l 

i=l 


with ?^e(n.*_ 1 ,r> i ) , n Q = A , h n+ 2 = B • from step 1, K( . . ) 
has a fixed sign. Hence det M can't be zero, which contradicts 
the previous conclusion that det M = 0 . This is true for 

ft _ ft 

any w +2 E ^ • Hence, either H^(w) = 0 , Vw £ I or H^(w) 
has at most n zeros (excluding the origin) in I . The former 
is impossible, because both H-^ and its argument identically 
zero in T , would require them to be identically zero V w . 


* 


AZ-23 


Step 3 . It is finally shown that 0(C) can be chosen so that 
H^(m) has n zeros (exluding the origin) in I so that 
in view of Step 2, it has precisely n . For C £ 1 j 

e(£) = nk( u ?-c 2 )<!>(?) (c 2 -a 2 ) (c 2 -b 2 ) , <{/(c) = I B.c 2Cn “ i:) 

1 1 0 1 

B q = 1 . (A4 ) 

B 

From (A3), h“(u). ) = / k n (w?-C 2 )tjKc)(C 2 -A 2 ) (C 2 -B 2 )dC . It 

2 A jj«i B 

is possible to choose iJj(C) so "that / (C 2 -A 2 ) (C 2 -B 2 ) C 2x ^(C )dC 

A 

for x a 0,...,n-l . For each x , this definite integral gives 
a linear equation in the coefficients B^ of ij^(c) , resulting 
in n linear simultaneous equations in the unknown Bj_ . It has 
been shown that this set has a unique solution [12, p.236], 
with n simple zeros of ij>(C) in (A,B) . 

For the present purpose, it is imperative to know the 
relation between the sign of H-^(w) in R k+1 and its sign in 
^k+2 *> see Figure A2), because in view of the previous four 
combinations shown possible for H-^, this relation must be 
flexible. By making use of the known result for the Tchebycheff 
system on [A,B] [11, p. 11] it can be shown that Sign(H^ in 
= (-1) Sign(H-^ in • Hence, to secure the 

desired flexibility it may be necessary to insert an extra 
zero say at the origin, or at A or B + , where it does no harm 
Of course, any desired sign of in is possible by 

choosing the sign of k in (A3). Also the scale factor of H 
is controlled by the magnitude of k , in order that L^H does 
not cross into B°(to) , V oi . 
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It is necessary to prove that l+LH has no rhp zeros. 

cl 

Obviously, 1+L has no rhp zeros, or L would not be 
an acceptable candidate for L ^ . The transition from L a 
to L H can easily be made continuous with L H £ B°(w) over 

a. a. 

the entire transition range, simply by means of the gain factor 

k in (A4). The zeros h l+L a H are continuous functions of 

k in H . Hence if these zeros are in the fthp at k = 0 

and do not Ife on the imaginary axis during the transition (by 

having L H £ B°(m) , Vw) , then 1+L H also has no rhp 
a. a 

zeros. The theorem is therefore proven, providing L a (joO 
has no values such as C in Figure A1 , V w 6 [O, 00 ) . 


If L (jw) = C (Figure Al) at one or more to values, one 
a 

may instead specify H^(m) of (Al) to be zero in I(A,B) and 
choose it in I so that ^(w) may have the required zero 
crossings in I . THe Hilbert transform giving H 2 in terms of 


H, 


(corresponding to (A3)), has the same total positivity 


property as (A3). This is satisfactory, providing L a (jto) has 
no values such as D in Figure A-^ , which is, in this respect, 
the dual of point C . In the most general case where L a (jto) has 
both kinds of values over one or more intervals, proof of 
the theorem requires a more general Hilbert transform - one 
in which Hg is specified in alter’nate intervals-say in 
1 / 


A 


(0, w a ) s I = (m b ,to c ) etc, j-h -*- b - 


I-. = (to ,to.) etc. 

D c d 

in Figure Al, Hj = 
E , F type . 


H, 


in Ip, = (w a ,to^) , 
is made zero in intervals of C type 


in the D type and either in the A, 


0 


a a -zsr 


Again, there is only one segment I = (A,B) in one of 

these intervals in which nonzero H 2 (£) Cor H^(O) may 

be assigned. A certain pattern of zero crossings is required 

in I . Figure A3 shows two examples, with the desired zero 

crossings and sign of in I A , part of I c , and I £ the 

same for both examples. The zero crossings for both the H 2 

shown are really same, but the desired sign of H 2 is the 

opposite in the second in the interval I_ . It will be seen 

later, that it is therefore necessary in the second example to 

introduce two more zeros, one at w + , the second at uk- , 

a b 

i.e. n = 11,13 respectively for the first and second examples. 


It has been shown how to derive the Hilbert transform for 
such problems in general and the solution has been given for 
the case of two intervals [7* 5, p.319j. For the special case 
here, where the specified function is nonzero only over I = (A,B) 
the same technique gives 


B ?H 9 Cc)dC A B . . 

! - J -, = vM i 


H, (to) = vCw) j — ^ ~ 

1 A (a) -C)A(C) 


a 2 r 2 

A oj -C, 


(-1) 0.5irH^(m), w £ I A • , 

« 

(-1)^0. 5 HgCw), m £ Ig,Ig, 


AC?) = IwCtu^-e 2 ) l 1/2 , v(w) 

A 


i = 0 for to £ I ; 1 for co £ I,-, 

3 L 


j = 0 for w £ I 


1 for w £ Ij ; 


B » 
d 


, 2 2 , , 1/2 
IT C CO “W J I 

A 


x " a , b , ... » 


(A5 ) 


Note that A(£) > 0 , V? 6 I J and v(w) > 0 , oj £ I . , I 

with v(w) ■*- 0 as w -> ao , <u K , . . . . The problem is to 

ci D 


choose 0(?) of (A5) so as to obtain the desired zero crossings 
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and signs of H^Cw) in 1^, I c , ... and of H 2 (w) in 
I I 

X B ’ D 9 “ * 


The proof for this general case is amost identical to that 
for Figure A2 and equation (A3). The ordering of the w^, rij 
is the same. Step 1 is exactly the same. Step 2 is also identical 

but note that the zeros referred to are the internal ones 

UptUj,. .. in Figure A3, not to the zeros at 

the latter being obtained by means of v (w ) in (A5). Step 3 
is also the same, as are also the remarks re R^+i > ^k+2 

which now refer to the relative signs of H^(w) at w = A- , B+ . 
(In Figure A3, n+1 is even and the signs are the same, so no 
additional zeros are needed for this purpose). However, one 
new complication concernes the relative signs of H,(u> -) , 

J. cl 

H 0 (w + } : of H 0 (w.-) , H, (w,+) , etc. It is seen that because 

l a l D X D 

of the factors (-1) 1 , (-1)^ on the right side of (A5) , no 
extra zeros are needed if Sg H^w^-) = Sg , 

S g H 2 (ui b -) = -Sg H 1 Cw b +) , Sg H 1 (oi c -) = Sg H 2 (m c +) . Sg Hj (« d -) = 
-Sg H^(to d +) , etc. An extra zero is needed for each violation 
of the above. Thus, in Example 2 of Figure A3, there are 


violations at , o> b 


and at w. 


(or w b +) 


requiring a zero at u>^ + (or tu^-) 
, giving n = 13. 


This completes the proof of Theorem 1. It has been shown that 

if a loop function exists at all which satisfies , then the 

optimum, if it exists, lies on at each w . It is implied, 

of course, that the B are bounded which explains the need 

u thihsrtidS 

for the supplementary ^boundary discussed in II 2. 
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APPENDIX 2 

EXISTENCE OF L . 

opt 

To prove existence, use is made of the analog of the 
Weierstrass theorem which states [13] that if JCL) is a 
continuous functional defined on the normal and compact 
family (L(s)} , then j CL) = min has a solution within the 
family. Here JCL) is equivalent to lim | ( jw ) e L( jtu ) | 

If the nominal L has no rhp poles then the family (LCs)} 
consists of functions regular in the right half-plane, 
including the jw axis which satisfy the bounds on L(jw) at 

each a) , and with lim arg LCjut) = -0.5Tre . If all PCs) 

QJ-^OO 

over the range of uncertainty have rhp poles then the nominal 

L must also have rhp poles 4 If at the nominal parameter set^, 

these rhp poles are at 1/t^,.,. 1/T m then let L = L^Cs)/ 

(1-t-j^s) ... Cl-T^s) and consider optimisation on L^(s) with 

the family {L^Cs)} regular in the rhp etc, and lim arg 

lim arg L, (jw) -*■ 0.5irCm-e) [5, Section 7.14]. If L has no 

poles at the origin then it is obviously possible to bound 

|LCj(u)| , (or } L-^ C j to ) j ) by some very ]arge number M , if 

necessary . If L has poles at the origin these poles can 

be removed and optimization considered on {sL} (or {s L^} ). 

If in (1), ACoi) = BCoj) is wanted at a finite number of values 

to a 3 then any candidate L would have factors 

2 2 2 2 

(s "Ho ^ ) ( s ... in the denominator. Such poles would be 

removed in the same manner. In Appendices 1-2, L is therefore 
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assumed mp and finite on the joj axis. By the principle 
of the maximum modulus |L(s)| (or | ( s ) | , or jrsL(s)|etc) < M 

over the entire rhp including the boundary. By Montels ' 
theorem [13], the family is therefore normal and compact. The 
functional J(L) is clearly continuous, so the conditions of 
the analogue of the Weirstrass theorem are satisfied and an 
optimum L(s) (or or s L etc) exists, providing the family 

{L} is not empty. { L) is certainly not empty in the mp 
problem [1] but it might be empty in the nmp problem. However, 
if the nmp problem can be solved at all , then the optimum 
solution exists and lies on B w at all w . From the nature 
of B y in the transition range (Figure 2) , l. t is clearly 
not realizable precisely as a rational function but may be 
approximated as closely as desired by a rational function. 


APPENDIX 3 
L opt IS UNIQUE 

Uniqueness of L ^ is first proven for the mp problem 
under the assumption that the B^O') are single-valued functions 
of 0' , the angle of the line to B^ from the -1 point. 

Conditions (Al) on B and the definition of the transition 

(JL) 

boundary (p, CM, -90e) at the end of Section II are important 
and should be noted. Suppose there exist two functions L a , 
which satisfy all the necessary conditions for optimality. 

Let to be the frequency at which L 'reaches' CM and similarly 
define . Let & = smaller of the sets [0,w a ] , [0,^] , 
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and define 

L ( jto)-L. ( jw) 

N(jw) = , 0 = Arg N , a. = Arg (1+L. ) 

l + L a (j<o) Nil 

ip(ti)) = Arg[L a ( jw)-L b ( jw) ] . ^ 


All candidates for L ^ are constrained to have the same excess 
(e) of poles over zeros and the same number (n) of poles at 
the origin. In the following development n = 0 is used, but 
if n t 0 replace L_^ everywhere by s n L^. Hence <3^ 0 >=0^(0) giving 
L^CO) - Lj^ CO) because both are on the unique B o (0) . 

Uniqueness is first proven on the condition L a , have 

no PM (phase-magnitude) crossing, defined as follows: 

L , L, have a PM crossing at u> iff 

^ ^ I_i ^ Li 

L a C ^ a 3 x ) = L b ( ^“x ) and s S n Ar S (rr) = “ Sgn Arg (it) 

^ b / to - v b'w + 

x x 

( A7) 

i.e. there is a zero crossing of Arg (L_/L. ) at w 

a. D X 

To simplify the proof, L = L, tangentially, i.e. with 
no zero crossing of Arg (L a /L^) also does not exist at any 
w 6 J) , although they can be arbitrarily close. Another assumption, 
later dropped, is that after L reaches the CM portion of the 
transition boundary ( Figures A4 , 5 ) , it remains on it until the 

O 

final portion Arg L = -90e (EDO in Figure A5a) . Suppose 
a (0 + ) > a, (0 + ) , then it is seen that a (w) > a, (to) and 

cl D a. D 

Qj^Cw) > 0 ,V m £ fi , as follows. The first is true because 
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\ 


5^(0' ) is single-valued, so a a = cx^ at to^ would require 

L a = at to^ , disallowed. To prove the latter, take 

as fixed with a (u> ) > a, (to ) . Inf 0. r (w ) , as 1 1+L, (jto )| 

a o DO No 7 D O 

ranges from zero to °° = Inf {Arg [(1+L )-(l + L K )] - a } = 0 , 

a d s. 

at |l+L a | = 0 which is an impossible value (the loop is then 
unstable). Hence 0 N > 0 , V to £ S' 2 . 

Consider (A) the case CM > 1 and ok < to (Case Al) at 

b a 

which, from the previous, ot a Cto ^ with ^(w^) of (A6) > it 

(Figure A4) . Suppose in Figure A4a L reaches CM 'behind' 

J a 

L. (Case Ala) meaning by this that Arg L, (jto ) < Arg L (jto ) . 

■D d a 3 . a 

Then L is 'behind' L K in this sense V io > to , because 
a d a 

PM crossings and departures from CM are not allowed. It is 

readily seen that as 0 )-*-°°. If reaches CM 'ahead* 

of L, (Case Alb) , i.e. Arg L K (jto) > Arg L (jto ) , the" it is 

seen in Figure A4b that i/i ( to = ) > 2 tt , L stays ahead for all 

to > to , and 0 N it as to -»■ » . The third possibility is 

(Case A2a) to^ > m a with L & reaching CM ahead of , 

giving 0 ^ ->• -it as u + « . The fourth (A2b) is > to and 

ahead of L a on CM , giving 0^ ■+ 0 as to + 00 . There 

are similarly four possibilities withing the alternative CM < 1 , 

of which 6^ -*• -tt or zero, as to ->■ , in three of them. 

But in only one case (to, < to and L. ahead of L on CM < 1) 

J b a b a 

does 6^j -+ -0.5ir e as to 00 . This is the only case in which 
the behavior of 0^ at infinity is consistent with the fact 
that from (A6), |N(jto)|- Oto -e inasmuch as and is 
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each Os~ e , and therefore if NCs) is mp and has no right 
half-plane poles, 6^(°°) = “O.STre . N(s) is obviously stable 
from its definition in (A3). It is conceivable that N is nmp , 
but if so 0 (°0 must be even less than -O.Sue , because 
right half-plane zeros effectively contribute [5, p. 332] phase lag 
at oo . It is therefore necessary to prove that this case is 
impossible, as follows. 

In Figure A5a at < w a > (ju^) is say, at point 

B on CM and at A must be to the left of O' X , 

because g q and a > in J2 . Both A and B must be on 

B (0 1 ) which is single-valued, and hence must be as shown. 
m b 

Since Lj, is 'ahead' of on CM , Arg < Arg , 

both on CM . The transition to the latter from the situation 

in Figure A5a with L (jw) in the interior of CM for 

to Cw, ,w ) , necessitates a (w ) = a. (ui ) for some w C (oj. ,to ) 
b a ao bo o £> a 

as shown in Figure A5b . But B must cross CM below B ' 

% 

(otherwise ^(jw ) would be in a forbidden region) and pass 

through A r . bote that for all u> e (o) K > w ) , B must 

Daw 

cross CM below L, (joj) , and in addition the B must be 

b id 

continuous (from Al) . These, two conditions, together with 

B of Figure A5a, compel B to have the shape shown in 
b a 

Figure A5b, violating its single-valuedness . Thus, this case 
is impossible. The conclusion is that if no PM crosssings 
are allowed of and each stays on CM once it 

reaches CM then there cannot .exist two such functions 
satisfying the optimality requirements . 
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Suppose that departures from CM are now permitted 
(but obviously only in the second quadrant). It is clear from 
the preceding that for CM < 1 the starting point is at the 
smaller of io_ , w. with a > a. and > a >0 at this 

a. D ci D a 

point. In order that 0 w C ^ > be consistent with |N(«0| , it 

is essential that (L -L.) execute a clockwise rotation while 

a £> 

satisfying all the constraints. The closest to such a rotation 

previously was in Figures A5a, b, invalid because point A' 

had to be on B in the interior of CM (otherwise, since 

% 

no departures from CM or PM crossings were allowed, L 

a 

always remained 'ahead' of , giving ii;(°°) = - it for 

e = 4- and in general*' -0 . 5ire +' ir ) . But now A^ can be on 

CM, but the resulting B must be as shown in Figure A6 (it would 

9 / 

be possible but more difficult to draw B if B' was on CM). 

r to 

o 

The desired rotation of (L^ - L^) requires that at some , 

L (jojj) depart from CM so that B must pass through 
a G Cd ^ 

L a ( ) j as shown in Figure A6 . Also, in order that the 

achievment at w not be lost, L, (j<o.) at B" must be in 

o b J d 

the shaded region and must of course not be in B^ , necessitating 
B as shown (it is easier to draw the latter if B" is CM ) . 


“d 

It is clearly impossible to continuously distort B into B 

w o 

. without crossing any continuous curve drawn between B' and 

o 


B" , thus forcing L b (jto) £ for some w £ (to o ,w^) . Hence, 

the desired rotation of (L^-L^) not ach i evable > even if 

CM departure are permitted. 


i 

i 

i 
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The same result is obtained in the case CM > 1 , where 
the best chance to achieve the desired rotation of L -L, is 

3 D 

when giving tt > tjj > 0 . 5 tt as in Figure A7. The 

zero crossing of from positive to negative is achievable 

if L is on CM - at say w . But the desired rotation 

requires that depart from CM while remains 

negative. must therefore advance clockwise, overtaking 

L^ on CM and forcing it off. However, in the process 

is forced to become positive again, in view cf the 

definition of p in 112. The desired rotation of L -L. is 

0) a b 

not achievable. Hence, L Q is unique if no PM crossings 
are allowed, even if CM departures are permitted. 

Let PM crossings of now be permitted. If there is 

a single one at to ,, £ 12 , the the roles of L , L. are 
° x a b 

simply interchanged in the proof above. This is so for any odd 
number of PM crossings in J2 , while the proof is unaltered if 
there are an even number of PM crossings in 12 . In the interval 
between w a , to ^ (or w ^ there can be no PM crossings 

because, by definition, one L is on CM while the second is' 
in the interior of CM . There can, however, be PM crossings in 
the interval (w,., 1 ®) where w is the larger of w , w, . In 
view of the previous argument , PM crossings in £2 need not 
be considered, only PM crossings in (w c ,°°) . At any such 
crossing | N| = 0 , 2,n|N| = -«> , which is accompanied by an 

abrupt positive change [5, p. 336] in 0^ . Thus PM crossings 
in (w c ,°°) increase the value of 0^ at Infinity, increasing 
even further the inconsistency between |N| and 0^ as w * 00 . 
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Proof of uniqueness is essentially the same in the nmp 
system. From Equation (3), the nominal loop transmission L 

no 

L no = A q (s) Lf mo C s ) is nmp because L mQ is mp while A q is 

all-pass nmp. The boundaries B n of L have the same 

w no 

properties as of L mQ in the purely mp problem 

A 

previously studied, i.e. 3 w" , 3 ,V m > oi" they do not 

enclose the origin and lie entirely in the second and third 

quadrants. Also they satisfy (Al) . The constraint that B n C Q ' ) 

to 

are single-valued in 9 ' is also meanwhile retained. N(joj) of 
(A 6 ) is defined exactly as before with nmp L , L, two 

cl D 

candidates for the optimum L , so that each must be of the 

no 

form L A (s), L mp . Thus N(s) has A (s) as a factor, 
mo m o 

The only difference in the proof is in the relation between 
jN(jw) | and e^Cw) as w -*■ . Since |n| is 0s~ e 

Bj^( 00 ) = -0.57Te -mt , where m is the . umber of zeros of A(s) . 
Hence the inconsistency between iN(°°)| and which was 

used to prove uniqueness, is magnified in the nmp system. 
Therefore in the nmp case, too, there is only one L function 
which satisfies the optimality conditions. 

Multivalued Boundaries B ( 0 ' ) 

to 

The final step is to remove the restriction that B (e 1 ) is 

tu 

single-valued in 0 ' . There is no difference between the mp 

and nmp cases here so the mp case and notation is used in 
the following. It is worth noting that this restriction of 
single-valuedness is actually a very reasonable one in practice. 
It is readily seen that to obtain multivalued B ( 0 1 ) it is 

M 

necessary to have plant templates (ranges of V in Equation (3)) 
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which can have at least four values for a given 0' , and 
system specifications A Jin J T j > A £n|L| and |~~| >:> ^ s 
thereby permitting B to closely follow the boundary of the 
plant template. The latter might be in special cases permitted 
at relatively large w , at which however, the plant templates 
tend to be vertical lines in the Nichols chart. Hence, the 
combination required for > multivalued B O') is highly unusual. 
Nevertheless this possibility is considered, for the sake of 
generality. 


Suppose that some B^ are multivalued in 0' , as in Figure 
A8 in the sector XY . Let a single-valued approximation B^ 
be made of B , differing from B only for some part of this , 
sector, e.g. may be taken as . . . QEFGH^JA^ . . . in Figure A8 . 

Let such closed portions of B ^ at any u> be denoted by A. 


U) 


i.e. in Figure A8 , A consists of GH-^J . Let A y be chosen 
so that the resulting satisfy (Al), This is certainly 

possible if the multivalued satisfy (Al). It has been 

shown that the optimum solution L ^ to the problem 

exists and is unique, if any solution at all exists. If this 
L ^ does not lie on A^ at any w (e.g. it is at A^ in 
Figure A8) , then it is also the unique optimum solution to the 
original multivalued problem. This is seen as follows. Suppose 
that B 0 consisting of . . .E'F ' IJ . . . had been used, instead of 
B , (and similarly for all other w in such a manner that 
conditions (Al) are preserved). Clearly the resulting ^optl 
would be precisely the same as before, as otherwise the 
optimum for the single-valued case would not be unique. In this 
manner every part of the multivalued boundaries in the XO'Y 


a z-3L 


sector can be given an opportunity to be part of a B . 

3 X 

without affecting the original unique L ^ • 

On the other hand, suppose that on some B , , L 

c aim ’ optl 

lies in the interior of EFG-type segment e.g. at A 2 in 

Figure A8 , but there is no w for which L ^ lies on the 

GH^J-type segment of B a i (0 • Then is again the unique 

optimum for the original multivalued problem, seen as follows. 

Let . . . QE 'EG 1 HGH^JAp . . (and corresponding ones at other w , 

in such a manner that conditions (Al) are satisfied) , be used in 

B a 2 • Then as before, ^optl would still be the unique 

optimum, the essential point being that a part of the boundary 

containing A 2 as an interior point, is common to B^ and B^ 2 * 

In this manner, all parts of EFGHIJ (except for the isolated 
1 1 ! 

points A 2 , A 2 ) have their chance to be parts of the boundary, 

without affecting the result. An infinitesimal modification 

of at A 2 would result in an infinitesimal change of L 0 ptl 

11 r 

and permit A 2 (or A 2 > to have its chance etc. 

Clearly, one need reconsider the problem only if for some id 
values L -ti lies on A^ , e.g. the closed GH^J -type portion 


, . As LCjto) and the A are continuous in w , and 

alto Co 


of B 

B , B are single-valued, L 


„ must be in a over some 

optl 0) 


finite interval [m , w „] • It is necessary to modify the B 

x y aw 

until there is no w for which L . lied on A * It will be 

opt w 

clear from L ^ how to proceed. For example, suppose ^ 0 ptl 
is in the A-^J-type protion in Figure A8, for w immediately 
less than oj x and more than u>y , then clearly the true L ^ 
wants to be in the JIH-type portion and return to the JA^typ'e. 
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The B should be modified accordingly (while preserving 

conditions (Al) etc) to permit iu to do so. If L traverses 

optl 

A , e.g. is in a JAn-type segment for cj < w and in a GF-type 
segment , then L opt wants to go right through the 

multivalued portion, and the B aw must be modified to permit 
it to do so. In this manner an L> 0 p t can be conceptually 
obtained which lies on no A^ except possibly at its boundary 
points. From the previous discussion, such an L .j. is unique. 

APPENDIX 4 

ON THE EXISTENCE OF A SOLUTION TO THE NMP PROBLEM 


As noted previously, there may be no solution to the nmp 
problem. It is convenient here to work with the mp L m0 of 
Equation (3) and the shifted (see Example in II and Figure 2 

Following the discussion there, it is shown here, that if an L 
exists which lies on the B , it must have turned the corner at 
V , before Arg L is zero. 

Form Equation (A3), it is easily found that 


j [R(w 2 )-R(m 1 )] = 


Co 


= / xce)xu>dc 

o 


5 


~{ nRCw 2 )-R(0)3-CR(w 1 )-R<0) ]> 

2 2 

£ (w, ) 

xco = — 

,2 2 * , 2 
(w^-'c )(a> 2 -? ) 


( A8 ) 


where R = Jin j L | , X = Arg L . It is assumed that the high- 
frequency boundaries are as shown in Figure 2, i.e. have 

the same shape but are shifted by Arg A(jw). Suppose that at 


hZ-3% 

> 

/ O 

3 ^ 1 ) - 0 j but |LCjaj 1 )| has not managed to decrease 
to the corner point V which would have permitted it to turn 
the corner and become negative in phase. Hence, for w > , 

A mu.it therefore be positi.ve. If a solution L exists, then 
at some w 2 > /l must cross from positive to negative values 
and remain negative, eventually reaching -0.5 ire . Hence, in 
(A8), , X(£) < 0 , £ E (w 2 ,«) = 14 , >0 in C(o^,w 2 ) = 13 , 

and may have both positive and negative values from zero to 
seme a (ID and must be negative from w to w, (12) . In 

(A8), A(D <0 in II, 12 ; >0 in 13, <0 in 14. Hence 
the net effect in each of 12, 13, 14 is positive, while in II 
there may be a net negative effect, which it will be soon shown, 
is overwhelmed by the others in any realistic problem. Hence 
| L( ) I > |L(jw^)| • However, | L( jto^) | > V , and since 

|L(j(i) 2 ) | > [L(jw-^)| there is no possibility whatsoever that 
the corner may be turned at w 2 > W 1 ’ un ^ ess 'the ^ha change 
in shape and move upwards in Figure 2 , as tu increases from 
to (U 2 . 

To prove that the effect of II even if negative, is over- 
whelmed by the others, consider R(w 1 )-R(to a ) using Equation A8 
with , w a replacing w 2 , respectively, with the same 

C(m 2 -wJ) 

X(?) but X replaced by A-,(£) = — s : ^ , whose , 

1 (tD-? 2 ) (coj-s 2 ) 

sign is opposite to that of A(s) in (w a ,w 2 ) , and the same 
in (0,m a ) , (u 2 ,«) . In (ti),,®) |A^(?)| < |A(e)| , while in 

(0,w a )j |\^(c)| > |A(?)| , which is the only interval in which 
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X(c)l(s) < 0 • Hence if this interval suffices to make 


/ X(c;)X( 5 )de; < 0 , it certainly makes / X(<;)X 1 C^)d^< 0 , 

0 0 1 

giving |L(jw a )[ < |L(jw^)| . This is inconceivable in any realistic 

feedback problem because it involves a region of positive phase angle 
in the low w range (unusual but not inconceivable) up to co , but 
in addition requires that at w significantly to a , the bounds on L 
require |L(jci),)| significantly^ (LCjm )| . This is therefore 

rejected. 

In practice in a specific problem, how may one determine apriori 
whether L can reach V before arg L is zero? A necessary 
condition is given here. Let L = L^M with Arg L^ = Arg ( jtu ^ ) = 0 
for oi > , 0^ = Arg M = 0 w £ , 6^ ^ 0 for w > , 

| MC 0 ) | = 1 . From equation A3, Sgn{ Jtn| H( j to ^ ) | -&n| M( 0 ) | } = 

Sgn / oo 0^(j;)y (?)dA with y(^) r [5 l"" 1 >0 in (0,ui^) , 


<0 in , so that |M(jw^)| > |K(0)j = 1 . Hence 

|L(jai^)| > |L^(jm^)[ . Therefore Ltjw^) certainly cannot reach 

the level V if cannot do so. One can therefore search 

for such a function. Define the optimum as that which 

satisfies B for w £ [0,w,] with minimum value of I L, (jw,)| , 

and the constraint ~ /l~] ( jo)-, ) = 0 , w 5 . Then it < n 

be proven by the same methods as in Appendix 1 that L t lies on 
for each w £ [0,w^]. Therefore in che appropriate w range, 
specify Arg L-^ to hhve the values dictated by the right hand 
vertical parts of ■ B h(j , and use a reasonable estimate of Arg L-j^ 
in the lower oi range. 
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It is fairly easy to answer the above question if a flexible 
computer program is available for determining L(jto) which 
lies on the . In a nmp problem for which no L exists which 

lies on the and goes to zero as w -»■ °° , such a program gives 

an L(jrn) whose phase becomes positive before L has decreased 
to V . The program must break down at some point thereafter, 

O 

if its angle at large w must be -90e 
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APPENDIX 5 

SOME TIME DOMAIN PROPERTIES OF NMP SYSTEMS 

This appendix presents some properties of nmp systems which 
are helpful in assigning realistic time-domain bounds and to their 
frequency-domain equivalent. An important effect of the all-pass 
factor A(s) = ( 1-x s ) / ( 1+xs) in T (s) = T (s) A(s) , T mp , is 
that of 'time-delay'. A(s) is a first approximation to the pure 

• ” 2 S T 

time delay e , reasonably accurate up to w « 2/x . This 

suggests and it is experimentally confirmed that if the bandwidth 
of the mp T^Cjtu) is < 2/x , then the nmp time response is 
closely that of the mp , delayed by 2x . The accuracy of this 
approximation increases with time and it is of course, not valid 
for t < 2/x , the step response being negative over approximately 
this interval. Similar time-delay estimates may be made for 
higher order all-pass functions, by finding the frequency range for 
which Arg A(jw) tox^ , x d the eq« ivalent time-delay. 

In the nmp system the step-response c (t) may begin 

with undershoot i.e. c (0+) < 0 . This is obvious from the 

n 

initial-value theorem, if A(s) has an odd number of rhp zeros. 

In fact, if the mp system has the typical continuously 

differentiable unit step response with c (0) = 0 , c (0+) > 0 , 

m m 

c («°) a positive constant, then the nmp system with step 

response C (s) = (l-x,s) ... (1-x s) C (s) with all x. positive 
n 1 x m i 

real, has at least x zero crossings and c n (0+) > 0 if x is 
even, < 0 if x is odd. This is proven by induction. 
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Let x = 1 , c m (0+) > 0 and then from the initial value theorem 

c n (0+) < 0 . For large enough t , c m (t) is negligibly small 

while c m (t) > 0 . Hence, at large enough t^ c n (t) > 0 and so must 

have at least one zero crossing. In the general case, write 

C n CS) = C1 -t 1 s) Cy(s) , so c n (t) = ' Suppose the 

proposition is true for x-1 , so that c^Ct) has a T least n-1 

zero crossings, while c^Ct) must have at least one zero crossing 

between any two of c^Ct) . Also Sgn c n (0+) = - Sgn c. : (0+) Ly 

the initial value theorem. Hence c (t) must have at '!.ea: ; t x 

n 

zero crossings. The proposition was proven true for x - 1 , so is 
true for all positive .integer x . 

The maximum magnitude of the negative undershoot in the nmp 
unit step response is of interest. For the case of a single rhp 
zero at 1/x , it was found experimentally that the following 
approximate relation gave fair results, M = O.lxw ^ , where w g 
is the -6db point of |T n (jm)| . Thus, to reduce M for a 
given x , one must decrease the system bandwith and thereby its 
speed of response. This empirical result can be theoretically 
justified by noting that in the small t interval where the peak 
undershoot occurs, the step response of the mp system .r T^ = 

T n = (1-sx) T m , can be well approximated by h m = At r , so 

. t 

h = h -xh c = f h d£-xh . At the instant t of peak 

nmmn^m^m u ^ 

undershoot c(t)=h(t ) = 0, giving t = n?' and 
n u n u ° u * 

M = -c (t ) = — — - — . In considering alternatives for h„ = At n 

n u , , ° m 

n+1 

it is reasonable to demand that at some t larger than t , 

o u 
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the mp 

Atr 1 . 

n+l 


t 

o 

step response / h d,$ have the same value K , i.e. 
0 m 

( \ n+ ^ 

K . Susbstituting in the above gives M = j n n 


5 


with n the only free variable, and M a minimum at n < 1 . 

This value of n is not feasible because n+l is the excess 
of poles over zeros of T^{s) and this excess must at least equal 
that of the plant. 
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£hp 

left half-plane 

Jtti 

linear time-invariant 

atv 

linear time-varying 

mp 

minimum-phase 

nmp 

nonminimum-phase 

rhp 

right half-plane 

















Example 1: 


Boundaries of L ( ju>) - Specifications I 

mo J 


HSSSESffliW 




Boundaries of L (jto) 

















• -Zero crossing 


Zero crossings of two complex building blocks. 
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APPENDIX 3 

Translation of Time-Domain Bounds on Unit Step Response 
Into Equivalent Frequency-Domain Bounds 

The results given here are largely empirical and were obtained by examining 
the works of many authors dealing with the problem and by considerable experi- 
mental work with many additional sample responses. Definition of the parameters 
of to and t-domain parameters are given in Figs. A3“la,b. 

Time delay = t ^ (Time for unit step response to reach O.l) 

Our empirical relation is 

to i t T ^ A ] (A3. ?) 

where A j is a function of the slope of | T ( j tu ) | for to > to j , the latter 
being the frequency at which the normalized frequency response is 0.1 (-20 db) 
in magnitude. The relation is shown in Fig. A3~2. This result is important in 
our design problem because it is essential that at high-frequencies the tolerance 
in | T ( jto) | be greater than the high-frequency plant uncertainty. This requires 
'widening' the nigh-f r rquency | T ( jiu) | bounds, and one should know that the 
price paid is in the time delay. 

Rise time, t “ t _ - t , (t Q = time for step response to reach 0.9) 
r o •_!_ • y 

Let be the average slope factor of J T ( j cu) ] for <o ^ > to > to_^ 

(defined by |T{jto^)| - “3 db, and to similarly defined). The results are 
given in Fig. A3~3- Among the samples used to obtain these results are those 


shown in Figs. A3-4 to 6 which vary very widely in shape. The agreement is 
wi thin ±10%. 

Relative Overshoot V 

A useful result is given by Papoulis [A1 ] for | T ( j ai) | having a single peak 
at Up and steady state response of unity, 

|T(j» )| 

1 + V < 1.18 ' | T (o) | ~ (A3. 2) 

Here we are also interested in the instant of peak overshoot denoted by t 

P 

The following empirical relation has been found by examining a large number of 
examples, ; or the case |T(jtii ) | 12 db, with accuracy of ±20%, 

(t “ t Q ) (w + (D ) fS K (A3. 3) 

P .9 p o s* p 

with K a function of |T(jw )| given in Fig. A3 _ 7- In Eq . (A3-3), <*> is 

p p o 

the value at which |t(jw) | = 1, normally u q = 0 in servo systems. The above 
Is valid for cases in which there is overshoot but no 'oscillation' (see below). 
Note that Eq. (A3. 3) Is used to find t , since t Q is known from the previous. 
Fig, A3“8 gives the peak response (1 + V) as a function of Aoj | T ( j to )|/u , where 

H r 

Au> = has been defined in Fig. A3“la. The accuracy of Fig. A3 - 8 is within 

±15%. 

Oscillations in step response 

A criterion for existence of 'oscillations' Is difficult to define. For 
example, we do not here regard (a) in Fig, A3 - 9 as oscillatory but do so for (b) . 
Roughly, there should be two discernible cycles with almost the same period. But 
of course we are here concerned with oscillations like (b) In Fig. A3“9, and not 


A3-3 


like (c) . The latter is associated with the peak at in Fig. A3 _ 10, the 

former with the peak at w^ . A great number of samples were taken with the 
result shown in Fig. A3~ll, where the 'peaking criterion 1 

A ^W. 

p.C. i—l |T(> )| (A3. k) 

P 

Auk = smaller of Auk, Auk , Aw. = w - w. defined in Fig. A3~la. From Fig. A3-11, 
J 1 Z 1 p I 

oscillations exist only if P.C. < ,07* 


Nonmin I mum-phase systems 

Appendix 5 of Appendix 2 presents some results for one or more right half- 
plane zeros in the system transfer function T(s) - Here we only present an 
empirical relation for the peak negative undershoot which occurs in the initial 
part of the unit (positive) step response, when the system has a single right 
half-plane zero at 1/a. The results are presented in Figs. A3^12a,b. Here 

refers to the -k db frequency for T(s) of the form T(s) - T^(s)- ^ ' 4." ' a ' s ' y t 
minimum-phase. The error A/MUS in the figures refers to the respective 
error obtained if the following relations are used 


MUS = (0.l)aw_ 6 
MUS = (.U)aw_ 3 

MUS is the maximum undershoot. 


(A3.5a,b) 
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